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Abstract 

The n-point correlation functions in single-field inflation obey a set of consistency conditions 
in the exact squeezed limit which are not present in multi-field models, and thus are powerful 
tools to distinguish between the two. However, these consistency conditions may be violated 
for a finite range of scales in single-field models, for example by departures from the Bunch- 
Davies state. These excited states may be the consequence of interactions during inflation, or 
may be a remnant of the era that preceded inflation. In this paper we analyze the bispectrum, 
and show that in the regime of theoretical control the resulting signal in the squeezed limit 
remains undetectably small in all known models which continuously excite the state. We 
also show that the signal remains undetectably small if the initial state is related to the 
Bunch-Davies state by a Bogoliubov transformation and the energy density in the state is 
small enough so that the usual slow-roll conditions are obeyed. Bogohubov states that lead 
to violations of the slow-roll conditions, as well as more general excited states, require more 
careful treatment and will be discussed in a separate publication. 
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1 Introduction 



Observations of the cosmic microwave background (CMB) and large scale structure (LSS) 
provide compelling evidence for nearly Gaussian, nearly scale-invariant primordial perturba- 
tions, just as one might have expected from a phase of inflation at very early times [1-4]. 
With the case for an early inflationary phase strengthening, it is timely to ask to what extent 
we will be able to answer more detailed questions about the very mechanism that produces 
it. One question that has gathered significant attention in recent years, because it is both 
interesting and seems answerable, is whether there was a single active light degree of freedom 
during this phase or many. 

An observable that is very sensitive to the spectrum of excitations during inflation is the 
squeezed limit of the three-point function of curvature perturbations. It has been shown to 
be of the form^ 

lim ia.a.aj = -p^{kL)Pdks) [k - 1) + o {ki/ki)] , (i.i) 

for models with a single light field [5-11], and more generally in dissipative models with 
multiple light fields but a single clock [12, 13]. 

The expression in Eq. (1.1), known as the 'consistency condition', is an expansion in 
powers of ki/ks which is valid to all orders in the slow-roll parameters. Its meaning is 
most transparent in the limit in which the wave-number vanishes, i.e. ki ^ 0. In this limit, 
Eq. (1.1) states (in co-moving coordinates) that the local physics is unaffected by the presence 
of long-wavelength modes. The fact that the leading corrections for small but finite values 
of kL scale as {ki/ksY is a manifestation of the equivalence principle. In other words, we 
measure deviations from free-fall through curvature effects. 

At first sight, the consistency condition appears very useful in distinguishing between 
different models of inflation. Any three-point function that does not behave as in Eq. (1.1) 
when ki/ks — )■ cannot derive from single-clock models. In practice, however, when we 
constrain models using data, we are forced to content ourselves with a finite range of momenta. 
For single-clock models, the behavior in Eq. (1.1) will eventually take over as ki/ks — )■ 0, 
but it may be violated for the finite k^/ks accessible to a given experiment. In order to rule 
out single-clock inflation experimentally, we must thus establish whether such models can 
produce a detectable non-trivial scaling for the three-point function for the observable range 
of momenta. Put differently, one needs to understand for which values of k^/ks single-clock 
models allow for deviations from Eq. (1.1). 

Scenarios that lead to a non-trivial behavior in the squeezed limit must involve either 
non-attractor dynamics [14, 15], or departures from the Bunch-Davies state. We will focus 
on the latter. There are two possible roads to excited states during inflation: either the state 
becomes excited by violations of adiabaticity in time-dependent processes during inflation, or 

^As usual ks and stand for the short and long wavelength modes, is the power spectrum without 
the i5-function, and (■••)' means we dropped the overall (27r)'^(5'^(fci + fc2 + ^3 + • • ■ )• 
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one assumes that inflation started in an excited state. If we assume tliat tlie excited initial 
state was generated by some unspecified but causal pre-inflationary dynamics, both scenarios 
lead to non-trivial correlations between long and short wavelength modes when all the modes 
are inside the horizon. In this regime, one can use fiat-space intuition to explain why single- 
clock models require an excited state to generate a non-trivial squeezed limit. Intuitively, since 
one can show that correlations in the Bunch-Davies state from early times are exponentially 
suppressed, excited states provide the necessary energy to produce non-trivial correlations 
that survive at late times. We will make this argument more precise throughout the paper. 

Whether obtained dynamically or postulated as an initial condition, the additional en- 
ergy density carried by the excited state will not only be responsible for the correlations, but 
will also affect the overall dynamics of the system. If it becomes too large, it may lead to a 
breakdown of the assumptions entering in our calculations, such as a meaningful separation 
into background and perturbations, or disrupt the inflationary dynamics altogether. Fur- 
thermore, the non-adiabatic time-dependent processes themselves must be sufficiently weakly 
coupled (since our universe is almost Gaussian). Otherwise very httle can be said about 
the resulting evolution of the curvature perturbations, and in particular whether there is a 
non-trivial scaling for the three-point function. 

In this paper we will take these considerations seriously and show that both for models in 
which the state is excited by non-adiabatic evolution, and for models with Bogoliubov initial 
states with significant occupation numbers but with energy densities small enough not to 
spoil the slow-roll conditions, non-trivial correlations can only exist between modes that are 
separated in their wave numbers by less than the fourth root of the observed power spectrum, 
or approximately two orders of magnitude, 

^ < A| ~ lO-l (1.2) 

We will discuss the consequences of this bound in detail. In summary, for models in which 
the state was excited by non-adiabatic evolution, it implies that they cannot produce a signal 
in current or planned experiments that are sensitive predominantly to the squeezed limit of 
the bispectrum. For the Bogoliubov states with large occupation numbers the bound (1.2) 
implies that they cannot lead to a power spectrum that is scale-invariant over more than 
two orders of magnitude. However, from measurements of the power spectrum from the 
microwave background, galaxy surveys, and the Lya forest, we know that the primordial 
power spectrum is nearly scale invariant over at least three orders of magnitude [16], thus 
ruling out Bogoliubov initial states with large occupation numbers. 

For Bogoliubov initial states with small occupation numbers and energy densities below 
the kinetic energy density of the inflaton, there are no constraints from the power spectrum, 
and the bound (1.2) is also weakened. We can thus not rule them out immediately. However, 
we will provide signal-to-noise estimates for Planck and Euclid showing that these models do 
not lead to observational signatures. 
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Once we allow the energy density stored in the excited state to exceed the kinetic energy 
of the inflaton, there are several things that require additional care. First, the energy density 
in the excited state is now large enough to affect the background evolution. In particular, it 
determines the evolution of H and since the energy density in the excitations redshifts like 
radiation, we now have a phase with 6 ~ —2. Second, on small scales k/a> 0^/^ the so-called 
(^-gauge is no longer well-defined. This is not too alarming, and of course also happens in the 
present-day universe, but requires some additional care. For this reason, we will limit ourselves 
to states for which the energy density stored in the excited state is negligible compared to 
the kinetic energy of the inflaton in the present paper and postpone a careful study of the 
regime in which this condition is not met to [22]. It should be noted that, strictly speaking, 
previous studies of excited initial states are also limited to this regime. 

States that are related to the Bunch-Davies state by a Bogoliubov transformation are 
superpositions of states with an even numbers of particles. As a consequence, the three-point 
correlations have to be generated through interactions, which are small in the case of slow- 
roll inflation. General initial states may have non-trivial three-point correlations even in the 
absence of interactions, and we will also consider them in [22]. 

In section 2, we illustrate both how the consistency relation can be violated for a finite 
range of momenta and what prevents us from making this range arbitrarily large. We use the 
lessons from this section to derive the constraint (1.2) for models in which the state is excited 
through non-adiabatic evolution as well as for models with Bogoliubov initial states with 
significant occupation numbers but with energy densities below the kinetic energy density 
of the inflaton in section 3. In Section 4 we generalize our results using the effective field 
theory (EFT) of inflation [17, 18]. However, the essence of the more general argument will 
be captured in sections 2 and 3. The more phenomenologically oriented reader, if satisfied 
with the exposition and results in these sections, may thus proceed directly to section 5 for 
an account of observational perspectives. We conclude in section 6. 

2 Origin of a non-trivial squeezed limit 

To understand what is required for a non-trivial behavior in the squeezed limit, let us first give 
a more intuitive argument for the validity of the consistency relation in the exact squeezed 
limit. Once the attractor has been reached, in single clock inflation modes freeze out after 
they exit the horizon. A frozen mode that is larger than the horizon is indistinguishable from 
a rescaling of the background and is therefore not measurable inside the horizon. As a conse- 
quence, the power spectrum of the short modes will not be modified at horizon crossing in the 
presence of the long wavelength mode. Therefore, if all correlation functions are determined 
when the shortest wavelength mode crosses the horizon, the consistency consistency will be 
obeyed. 

To evade the consistency condition in single-clock inflation, a correlation between long 
and short modes must have been generated before the long wavelength mode crossed the 
horizon. In the exact squeezed limit, i.e. for /c/, ~ this cannot happen because the long 



4 



wavelength mode is always outside the horizon. This is the hmit in which the consistency 
condition apphes. However, for finite values of momenta it need not be the case. Before 
horizon crossing, the long mode is distinguishable from the background and can induce cor- 
relations with the short wavelength modes. This correlation acts as an initial condition for 
modes observed a late times, and therefore it is observable. 

It is well-known [5, 7] that the three-point function computed in the Bunch-Davies state 
does not depart from the scaling in Eq. (1.1). What prevents such correlations from occurring? 

The answer is that, in the Bunch-Davies state, the contributions to the bispectrum from 
modes that are deep inside the horizon {k/a ^ H) are exponentially suppressed. We can 
understand the origin of this suppression from fiat space intuition since the modes inside 
the horizon are essentially insensitive to the curvature. For illustrative purposes, let us then 
consider in-in correlation functions in a time independent system. For example, if we have 
a scalar field with an interaction Hamiltonian ifint = l/^V^^, the equal time correlation 
three-point function is given by 

in^n^nM' = -^ra^ f rftV^^'^^^-^-^'^^)^*'-*) - c.c. (2.1) 

J — oo(l — ie) 
1 

^Akik2h{ki + k2 + h) ' 

where ki = \ki\ is the energy of the mode with momentum fcj. (Note the parallelism with the 
computation in an expanding universe, in particular the factor of l/fc*, with kt = ki + k2 + k^. 
The main difference, of course, is the late time behavior of the correlators in de Sitter, where 
curvature effects start playing a role.) Although it appears as if the correlation function 
receives contributions from all times, the contribution from early times is exponentially sup- 
pressed due to the rapid oscillations. 

One can intuitively understand the origin of the correlations, and the fact that contribu- 
tions from early times are exponentially suppressed, as follows. The computation that leads 
to the result in (2.1) represents the production of three virtual quanta from the vacuum at 
some point in spacetime that are subsequently "measured" by the field(s) fj^_(t)- We describe 
this process diagramatically in the left panel of Fig. 1. Since all three quanta originate from 
the same event, they are correlated. Because they carry energy, energy conservation prevents 
them from going on-shell, implying that they could not have been produced long before they 
are measured at time t. However, energy conservation may be violated quantum mechanically 
for a period of time At ~ 1/AE = l/kf. We conclude that the quanta must have been created 
less than At prior to t, explaining why the contributions to the integral from early times are 
exponentially small. This intuition implies that also in the Bunch-Davies state correlations 
are generated near t' ~ t — At. When translated to (, this implies that the correlations are 
generated shortly before the short mode freezes. Since the long mode is already frozen at this 
time, this then ultimately turns into the consistency condition given in Eq. (1.1). 

The reader might complain that the case of the three-point function is not very generic. 
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Figure 1: Sample processes responsible for correlations in the three- and four-point 
function given in equations (2.1) and (2.2) 



because there is only one possible way these quanta can be created. To address this, let 
us consider the four-point function {ip^^ip^^ipj:^^ip^^{t))' as a less trivial example. One of the 
processes that is responsible for these correlations is represented in the right panel of Fig. 1. 
Three virtual quanta are created from the vacuum with momentum (fci, /c2, q = — fci — at 
time tl . Two of them are measured by two of the fields at time t, while the third decays at a 
time t2 — tl + Ati (with Ati ~ l/{ki + k2 + q) in agreement with the uncertainty relation) 
into two quanta with momentum {ks, k^) which are subsequently measured also at time t 
by two other fields. Since the presence of the four quanta violates energy conservation by 
AE = ki + k2 + k^ + k4^ = kt, the instant ^2 should be within At = 1/kt of the instant t at 
which they are measured. This process thus results in the contribution to four-point function 
given in the first line below 



lQkik2k^k/i q 



{ki + k2 + q) k 



11 1 

+ V-, ; ^T- + 



(fcs + k4 + q) kt {ki + k2 + q) (fcs + k4 + q) 



(2.2) 
+ permutations . 



The second line arises from similar processes. We see that also for higher n-point functions 
correlations generated at early times are exponentially small, simply because quanta had to 
originate from the vacuum state and energy conservation prevents them from going on-shell. 
The more energetic the quanta the closer to the measurement at time t the correlations are 
produced. This also implies that these correlation functions do not exhibit non-trivial be- 
havior in the limit as one of the external momenta becomes soft. Note, however, that there 
is a pole in this correlation function as g — > in the internal line (because the particles are 
massless). This pole ultimately leads to the consistency condition for the four-point function 
in the collinear limit. (See [24-29] for more details on soft and collinear limits.) 

As we emphasized, the consistency relation must hold once the long wavelength mode 
is outside the horizon. Therefore any violation must arise while all the modes are inside 
the horizon and our fiat space intuition (and notion of energy conservation and uncertainty 



6 



principle) applies. In what follows we will then use this terminology loosely while referring 
to the possible realizations of a non-trivial scaling in the squeezed limit. There were two 
important assumptions that lead to the exponential suppression: a) the system is time trans- 
lation invariant, and b) we are computing vacuum correlation functions. If either one of these 
assumptions is violated, then the correlations are not necessarily suppressed. This is because 
energy conservation does not prevent the configuration from being realized. In the context 
of inflation, we can easily break either of these assumptions. In our flat space example, we 
might imagine an interactions of the sort ii{t)(fi^, with ii{t) — cos{ki,t) arising the time evolv- 
ing background field (f){t). Now the contributions to the integral are not suppressed when 
kt = ki,. This is consistent with energy conservation, as we have extracted energy k^, from the 
background (pit) to produce real quanta. In this simple example, the reader may recognize 
the possibility of resonances and oscillations in the power spectrum. 

The second assumption, that we are in the vacuum, may also be violated during or 
prior to inflation. Although the energy density will redshift away during expansion, at the 
onset of inflation it is plausible that the fluctuations were in an excited (i.e. flnite energy) 
state. Energy conservation no longer forbids contributions to the correlations from early times 
because we are starting in a flnite energy state, again, given by some scale k^,. In the flat 
space example, if the mode ki is excited, a negative frequency —ki can appear in the integral 
as e**^*, with k = —ki + k2 + ka. Now we get contributions at all times when k ^ 0, namely 
ki = k2 + k^. The reader may already recognize this as the contribution from excited states 
to the so-called folded shape. 

An important feature of both the exceptions we mentioned above, is that they are asso- 
ciated with some positive energy scale, i.e. /c*. For the case of departures from the vacuum 
state, there will be an energy density present from the beginning of inflation which will back 
react on the geometry, and therefore it must be bounded in order to retain the basic features 
of the inflationary dynamics. On the other hand, for the case of resonances, the scale k^, 
cannot be too large or else interactions induced by the non-adiabatic evolution may become 
strongly coupled, in which case very little can be said about the resulting evolution of the 
curvature perturbations, and in particular whether or not a non-trivial scaling develops. We 
then conclude that the necessary condition kt < k^ for a non-trivial scaling in the squeezed 
limit must itself be bounded. We will show how these considerations lead to an universal 
bound on A;^, hence to a lower bound on ki/ks- 

One might ask about the validity of the consistency relation when correlations are intro- 
duced as initial conditions before inflation. In this paper we will assume causal pre- inflationary 
dynamics implying that modes inside the horizon H{to) are uncorrelated with modes out- 
side the horizon at that time. This implies the validity of the consistency condition when 
kL/ks < a{to)H{to)/k^. 
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3 Bounds on squeezed limits in slow-roll inflation 



Let us now derive a bound on the range of wave numbers that can be correlated in a non-trivial 
way for inflationary models that (at least on average) satisfy the usual slow-roll conditions 

H H 

e = -— <1 and 5 = (3.1) 

2HH 

We will first discuss models in which the state is being excited continuously during inflation 
and then turn to slow-roll models in which the state was excited at some fixed early time 
by some unspecified dynamics (plausibly one which violates the slow-roll conditions). The 
discussion of states that are sufficiently excited so that at least for some period of time 6 ~ —2 
will be postponed to [22]. 

As we already explained, correlations between a long mode with wave number ki and two 
short modes with wave numbers ks must have been generated before (or around) ki ~ a{t)H 
when ks ^ a{t)H. We also know, from energy conservation, that at any given time only 
modes with wave numbers up to fc^ < fc^ are significantly excited and can contribute to non- 
trivial correlations at late times. Our goal will be to derive an upper bound on fc^, hence on 
ks/ki- In the case of non-adiabatic evolution, it will correspond to the validity of a weak 
coupling expansion for the perturbations. In the case of the excited initial state, it will follow 
from the validity of the slowly-varying background evolution conditions. In either case it 
turns out to be of the form 

kja<<p'/'. 



Therefore, if modes are significantly excited, non-trivial correlations can only arise between 
modes separated by 

''">H/^'/\ (3.2) 



kL 



We recognize the right hand side as the fourth root of the power spectrum^ 

^'(C.<_.^)' = ^? = ^ " 4.7 X 10-«, (3.3) 
which allows us to re-write our bound as 

1^ > 4 - Af ^ 10-. (3.4) 

The above expression for the power spectrum as written is only valid for the Bunch-Davies 
state, but we will show in what follows that in essence our simple line of arguments remains 
unchanged. 

As a result, COBE normalization puts an upper bound on the range of wave numbers 



^Here = 2tt^A^, where A| « 2.4 x lO^^ [1]. 
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ks/ki that can display non-trivial behavior in the squeezed limit. 



3.1 Non-adiabatic evolution 

Let us now consider the case in which non-trivial correlations are generated by interactions 
during inflation in more detail. 

In this case, the Hamiltonian typically has a significant time-dependence on a time scale 
At ^ uj^^ ^ H~^. As long as the natural frequency of the modes is larger than this frequency, 
k/a> w*, they evolve adiabatically and remain in the Bunch-Davies state. However, as they 
redshift their frequency eventually becomes comparable to and the state becomes excited, 
leading to non-trivial correlations. Concrete examples that we now discuss are resonant 
models [30-32] or dissipative models [12, 33-35]. 

3.1.1 Resonances 

As a first concrete example, let us consider a model of single-field infiation with the potential 
[31, 32] 

V = V^ + k^coJ'^\ , (3.5) 



where Vq is a potential that supports slow-roll infiation, and A, / are energy scales. As the 
infiaton slowly rolls over the oscillatory potential, the Friedmann equation implies that the 
background geometry will oscillate with frequency ^ 4'/f-'^ We assume that the modes 
start out in the Bunch-Davies state at early times when k/a ^ u^, ^ H. As the modes 
redshift their natural frequency becomes comparable to the frequency of the background and 
they eventually undergo parametric resonance and the state becomes excited. This intuition 
is in perfect agreement with the solution of the curvature perturbation 

Cdt) = Ckit)a{k) + Q{tWi-k). (3.6) 
Up to slow-roll suppressed corrections, the mode function after the resonance is given by 



k 



X) 



(3.7) 



where x = k/aH. As anticipated, the non-zero Bogoliubov coefficient /3, 



k 



TT 1 /o . , UJ* , , A^ 



/3k = 36oa/ 7t" ^^^e ' f with a = t7 and bo = , , (3.8) 
V 2 H K(0o)/ 

shows that the resonance has lead to an excited state. Here (pQ and (pk denote the value of 
the scalar field when the pivot scale fco and the mode with comoving momentum k exit the 
horizon respectively. Since \(3k\ <^ 1, the normalization of the power spectrum is thus still 



■^Note that for a slowly rolling (f> this is reminiscent of the time-dependent coupling we discussed in flat 
space, i.e. /x(t) ~ cos(w*i). 
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given by the slow-roll result 



tt4 

= ^ ~ 4.7 X 10-^ . (3.9) 



Formula (3.8) also reveals that the phase of Pk is oscillatory as a function of k (rather than 
constant, which is often assumed for excited initial states), which will be relevant to the 
observational consequences discussed in section 5. 

The bispectrum in the resonance model receives its largest contribution from the inter- 
action [30] 

Hint = - J d^xdta%t)e{t)5{t)C^C ■ (3-10) 

It gives rise to [32] 



(C.^.C,^.C,^3)' = [sin (alog(^^; 

Aji thn ri/Q L \ 



1 k 

a 



5^^cos(alog(^^±|f^)) +0(a-^)] , (3.11 



where the amplitude is given in terms of the underlying parameters of the model as 
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To see that there is indeed an intermediate regime in which the consistency condition is 
violated, let us first consider the behavior of the bispectrum for squeezed configurations in 
which ki <^ k2jk'j, satisfying a'^ < |^ ^ 1 with k^ = ki and ks = {k2 — ^3)/2. In this 
regime the term on the first line dominates over the term on the second line. Not also that 
in this regime we cannot Taylor expand the sine in (3.11) in a series in ki because > 1. 
As a result, we find B{kL, ks, ks) oc l/k1 which violates the consistency condition. Now 
consider configurations with |^ <^ a~^. For these modes the second term dominates and 
we can expand the oscillatory term. We find that in this regime the consistency condition 
holds [32]. By going to higher orders one finds that the leading corrections go as {ki/ksY as 
advertised [8]. 

The basic result here is not surprising, cu^ ~ 0// sets the energy scale of the oscillatory 
behavior which controls the largest values of ks/a{t) with which quanta can be produced. 
As a result, using k^ ^ a{t)H, we find that we need ^ < ^ if we want to violate the 
consistency condition. To make the regime in which the consistency relation is violated as 
large as possible, we must therefore increase u^, at fixed H. The value of is fixed by the 
slow-roll potential, and we adjust / to do this. However, theories with an interaction as 
in Eq. (3.5) are a valid effective description for scales below the cutoff Ajj = Airf. The 
energy scales we study should all be below this cut-off [45]. In particular, we must satisfy 
1^* < A[/, which implies that > and hence cj* < 0^/^. Once again, from the conditions 
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ks/a{t) < k^/a{t) ~ UJ^, and > a{t)H, we obtain 

|<^^:0^. (3.13, 

Intuitively, from our discussion in Section 2, we understand that to produce a squeezed shape 
we are trying to create modes with energy E ^ H, hj extracting it from the kinetic energy 
of By tuning the frequency in the cubic couphng, i.e. 6(t) ~ cos(a;*t), we can make this 
conversion possible for modes with higher momenta but we find that, for the validity of the 
effective description, we are ultimately limited by the kinetic energy density of the inflation, 
0^/^, which is the same scale that appears in the power spectrum. 

3.1.2 Dissipation 

Another example of non-adiabatic evolution arises when we allow to interact with addi- 
tional (dissipative) degrees of freedom [33-35]. In those scenarios the inflaton produces a 
steady burst of particles as it (slow) rolls during expansion. The dissipative effects allow for 
inflationary dynamics for a steep (er) potential. Interactions with these extra degrees of free- 
dom will then excite the inflaton significantly, from the Bunch-Davies state to a semi-classical 
state. (For example, this occurs in warm inflation for T ^ H.) Even though in general 
dissipative effects require multi-fields, one can show that the consistency condition is obeyed 
for a vast class of models, provided these additional degrees of freedom do not contribute to 
curvature perturbations [13]. 

The case of dissipation resembles that of excited initial states, although not as an initial 
condition but as a (recurrent) dynamical process, and one gets contributions to the three-point 
function that peak in the folded limit [12]. However, since dissipative effects are more efficient 
when the 'friction' coefficient 7 is much bigger than Hubble, the memory of the correlations 
is quickly lost and most contributions occur just before the modes freezeout. In practice this 
means there will not be a significant departure in the squeezed limit from what occurs in 
single-field inflation in a Bunch-Davies state. This was explicitly shown in [13]. 

3.2 Bogliubov initial states 

Let us now turn to excited initial states. As a cautionary remark, let us point out that for 
any model in which the inflationary phase lasted longer than the minimum number of e-folds 
required by the horizon problem, observationally we will not have access to the modes that 
were significantly excited at the beginning of inflation. So the discussion in this subsection is 
only relevant for models in which inflation lasted just long enough. Furthermore, for reasons 
we explain below let us restrict to a commonly studied but rather special class of excited 
states that can be obtained from the Bunch-Davies state by a Bogoliubov transformation. 
More general excited states will be discussed in [22]. 

In these states, it is convenient to write the curvature perturbation as 
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where 6t and 6g are creation and annihilation operators and 6^ annihilates the state. The 
Bogoliubov coefficients satisfy — = 1 and without loss of generality, we may take 

= ^1 + /3|e*^fe with real Up to slow-roll corrections, the mode function is given by 

a(r) = ^^(l + zA:r)e-^^ (3.15) 

As long as we limit ourselves to states with energies that are small compared to the kinetic 
energy of the inflaton, it is straightforward to compute correlation functions using these 
states [36-42] and the power spectrum is 

k'UM = (^c)' = ^ + + ^^\Pk\' + m'^os{e^)) ■ (3.16) 

The bispectrum is also readily calculated, though the final answer turns out to be some- 
what lengthy [37-42]. For simplicity, we only present expressions for the squeezed limit for 
large and small /3g as well as 6'^ ~ 0.'* For large the squeezed limit of the three-point 
function becomes 

where tq is the conformal time at which the state was excited. These expressions are useful to 
gain some intuition about the behavior of the bispectrum in the squeezed limit, but we will 
work with the full expression when we turn to forecasts for future experiments in section 5. 
We immediately see that we get an enhancement in the squeezed limit, provided ^2, /ca < k^,. 
Notice also that the result is slow-roll suppressed. This can be understood from the fact that 
these states are a superposition of states containing only an even number of quanta. As a 
consequence, there are no three-point correlations built into the state and slow-roll interactions 
are required to generate them. This will not be true for general excited initial states making 
the Bogoliubov states somewhat pathological, if the goal is to discuss the possibility of a 
non-trivial three-point function. We postpone the study of general initial states to [22]. 

3.2.1 Bound on ks/ki for Bogoliubov initial states 

Let us now derive the bound on ks/ki- From equation (3.16) we see that the requirement 
of scale invariance imposes /3g ~ /3g for (as well as a scale independent phase ~ 6^^ ) for 
all the modes we observe. Moreover, maintaining a finite energy density requires |/3gp oc A;~" 
with n > 4 when k > k-^. 

*A constant phase seems rather unphysical and does not arise in any microscopic model we know. However, 
at least for /3g 1 this is required by scale invariance of the power spectrum, and we will allow for this 
possibility. 
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As explained before, the constraint on comes from bounding the total energy density 
in the state 



Here Oq corresponds to the value of the scale factor at the time Tq when the modes are excited. 
If the energy density in this state is to leave the background evolution unaltered, it must be 
less than the kinetic energy of the inflaton. To see this explicitly, consider the equation for 
H in the presence of this energy density 

-M,^,ij^i0^ + i3i,|(^)^/3,J^ (3.20) 

So unless the second term on the RHS of this equation is subleading to the first, H is affected 
by the excited state, and to avoid this, we should impose 

(3^21) 

Using the power spectrum in (3.16), the above inequality becomes 
kt 67rV „ 2 / A/3 



1 + 2|/3,H 




p + m 


|4cos(^,^) 




2 



« rr.^o 12 ^c) — ^irn^ — — —■ (3.22 



The simple lesson here is that the scales that control the power spectrum, namely and 
Mp[|if| ~ (j)"^, are the same scales that appear in the back-reaction constraint on (> ks). 

Since ks < fc* and ki > aoH, we see that the constraint that the effect on the background 
evolution be negligible has enforced 

^^<102 for /3,^^»1, (3.23) 



L 



k^ 

ks < 10^ 



for /3g <1. (3.24) 



Large occupation numbers are thus immediately ruled out because this bound would imply a 
step in the power spectrum, which is not observed. 

For small occupation numbers the bound is weakened and occupation numbers around 
Pk ~ 10~^ certainly are compatible with the observed power spectrum. However, combining 
our bound (3.24) with equation (3.17), we see that the amplitude of the three-point function 

1/2 

in the squeezed limit decreases hke /3/ as we decrease /3r . One might thus suspect that the 
amplitude will be too small to be detectable. We will be more quantitative about this and 
discuss the observability in more detail in section 5, showing that this suspicion is correct. 
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3.2.2 Generalization to small speed of sound 

One can easily generalize our previous argument to the case of small speed of sound models. 
The main difference is the computation in Eq. (3.19) for the energy density, which is now 
given by 



(27r)3 

As in the case of slow-roll, this contribution must be negligible compared with M'^^H . As a 
result, we find 

<ImI\H\ . (3.26) 

The condition for freeze-out is modified due to the sound horizon, and therefore we now require 
ki > aQH/cg for the long wavelength mode to be inside the horizon when the correlations are 
created. Using the power spectrum, and following the same steps as before, we obtain 

' < J ^ ^37rMAg c^ (3.27) 



kt H^\P^,f - V C 

If we demand that there is no superluminal propagation, namely < 1, then the bound on 
ks/ki is at least as strong as before. 



4 Bounds on more general models: The EFT of inflation 

In this section we explore excited initial states and non-adiabatic evolution from the per- 
spective of the EFT of infiation, developed in [12, 17, 18, 43]. Here we will only cover the 
necessary steps to generalize the results from our previous section to include a wider class 
of models, for instance those with a non-trivial dispersion relation. We will elaborate on the 
more formal aspects of the EFT formalism elsewhere [22]. (See also [44].) 

The EFT of infiation is characterized by the spontaneous breaking of the time diffeo- 
morphisms of the FRW background geometry. As a result, one may write terms in the action 
that depend explicitly on time, provided the diffeomorphisms are realized non-linearly. This 
can be achieved by introducing the Stiickelberg field vr, and writing the action as a function 
of t -|- TT such that it is now invariant under: t — )■ t + ^°(a;, t), tt — t- vr — ^°(a;, t). The action can 
then be written as 

S = j (fx^ [ \MlR - F{t + tt) - c{t + ■K)d^{t + Tx)d^{t + vr) + ■ ■ ■ ] . (4.1) 

The functions F{t) and c(t) are constrained by demanding that the background geometry is 
an FRW universe with metric 

ds^ = -de + a{tf ( + r^dnA . (4.2) 
y 1 — nr^ J 
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It is straightforward to solve Einstein's equations. For situations when we only have one 
contribution to the energy budget of the universe one finds [17, 18] (setting k = 0) 



c{t) = -Ml,H (4.3) 
F{t)=Ml,(?,H\t) + H{t)) . 

In the presence of an excited state, the energy density stored in the excitations, i.e. p-,,, will also 
contribute to the energy budget and background equations, and may even fluctuate. In such 
a scenario one should distinguish between two type of perturbations, those associated with 
the 'adiabatic mode', which we denote as vf, and those which account for the perturbations 
of the clock that controls the end of inflation, which we will retain as vr. For the former (4.3) 
is still valid, whereas for the latter we have [12, 13] 

c(t) = -M^,H-lp^{t), (4.4) 
Fit) = Ml (?,H\t) + H{t) - |p.(t)) . (4.5) 

In this paper, we work in the regime p-j^ <^ Mpj|if | so that the two are essentially the same 
TT ~ TT. The regime pj^ ~ M'^^\H\ for which the distinction is important will be studied in [22]. 

4.1 Non-adiabatic evolution 

In section 3.1 we discussed resonant non-Gaussianity which, in the language of the EFT of 
Inflation, arises from the time dependence of H{t) ~ -f^o(^) + -f^osc(i) cos(a;^(t + tt)). The 
key feature that leads to a non-trivial squeezed limit was the rapid time dependence from 
Ui, ^ H . We can now generalize this analysis to any form of non-adiabatic evolution by 
considering generic time-dependent couplings. We will focus here on models that lead to 
non-trivial squeezed limits and study the constraint on the number of modes involved from 
requiring the EFT is weakly coupled. 

Consider adding a contribution to the Lagrangian for tt of the form 

5£int D MV(t + vr) g{d^{t + vr), ^^(t + tt), ■ ■ ■ ) , (4.6) 

where f{t) ~ u^,f(t), g{df^(t + Tr), 9^(t + 7r), ■ ■ ■ ) is polynomial in its arguments and M is some 
mass scale required by dimensional analysis. In the spirit of generality and due to the natural 
features of the models, here we will deal with an energy scale rather than a momentum 
scale fc*, as we did in section 3.1, to allow for generic forms of the dispersion relation, i.e. 

= Cp{uj^)h. 

We will be interested in the case where the correlation functions are nearly scale invariant, 
and therefore impose an approximate (discrete) shift symmetry: f(t + c) ~ f(t) [45]. As a 
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result, we expect that at a generic instant of time: 



£r/(t)~u;:/(t). (4.7) 

These criteria imply that the scenario we have in mind is effectively a generalization of reso- 
nant non-Gaussianity [32] (some such generalizations were discussed in [45, 46]). 

We will take the kinetic term for vr to be of the form £kintetic = A^tt^, such that the 
canonically normalized field is tTc = A^vr. (In slow-roll inflation, one has A'^ ^ cfp'.) Expanding 
Eq. (4.6) in powers of tt then we find 

5£iNT 3 M'ujy^f{t) g{d,{t + vr), d'it + vr), ■ ■ ■ ) = g{d,{t + tt), ^^(t + vr), ■ ■ ■ ) 

(4.8) 

In the limit of large n, the details of g are irrelevant for determining the strong coupling scale 
Astrong ~ A^/o;*. As in section 3.1, the contributions to the squeezed limit arise at energies of 
order w*. 

4.1.1 Generalized bound on ks/ki 

In order for the expansion in (4.8) to be under control, we must require 

« 1 . (4.9) 

If this inequality is violated, the effective theory is strongly coupled at the scale o;^ and 
therefore the bispectrum is incalculable. Therefore, to generalize our results from (3.1) we 
need to place an upper bound on A. If there is no such bound, one could make cj* arbitrarily 
large without violating (4.9). Following [47, 48], one can show that the null energy condition 
implies 

|if I > A^vrl^^,,, = A^vTeU^^, > w'vr'U^^. , (4.10) 

where in the last step we used (4.9) and vie = ^t^c- Our theory is nearly Gaussian with a 

3 _ _3 

kinetic term tt^, therefore vTc ~ /c 2 a; ^'^ = Cp at high energies, which leads to 

Ml\H\ > ujjc^'M . (4.11) 
Therefore, correlations between short and long momenta can only occur in the range 

k, - c,{co ^ CO.) H H cy\uj.)^ ' cyv)~ ' ^ ^ 

where we used (see appendix A) 

Ml\H\/H^<\Afc;\uj^H). (4.13) 
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In the final step, we have ignored the ratio of phase velocities, which is justified by the small 
fractional power and that the phase velocities grow with energy for simple dispersion relations. 



4.2 Bogoliubov initial states 

Let us now generalize our bounds derived in section 3.2. As before, we restrict ourselves 
to the case in which the energy density of the excited state does not affect the background 
evolution. From equation (4.4), we see that this implies 

p.(t)«|M2|if|. (4.14) 

4.2.1 Generalized bound on ks/ki 

Starting from (4.14), we can now follow the steps in 3.2 to write this inequality as 



a(ro)3 J (27r)3 ^'^(-o)^'^^' 



(4.15) 

The main difference with Eq. (3.19) is the scaling, because here we did not assume a dis- 
persion relation of the form w ~ fc. We can now show how our bound holds more generally. 
Specifically, we can always write (4.15) as 

H'cl(Kr„)\^i-\ < 12' — . (1-16) 
where Cp(A;, r) = '^j^j^ is the phase velocity. Moreover (see Appendix A) 



<i(^c)~ -H)[l + 2|/3,H^ + + cos(%)] . (4.17) 



This inequality implies that the right hand side of (4.16) is bounded by the power spectrum. 
Using once again ks < A;* we obtain, for p > 1 (and barring tuned cancellations in the 
power spectrum), 

^;|g^< 6.= (a?) (4.18) 

We can then use this inequality to generalize our bound for the ratio of scales: 

ks _ cM Uj{ks) ^ /C^X 1/4 / -1/2 ^ ^^^^g^ 



ki Cp{ks) H Vcp(/c5)7 V 

Hence, even in the presence of a non-trivial dispersion relation, the bound remains essentially 
unchanged. 

The reader may have already noticed that the bound in (3.27) was somewhat stronger 
than the one above. In other words, (4.19) does not reproduce (3.27) when Cp = Cg. The 
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reason is that here we used inequahty in (4.17), which is vahd regardless of the form of 
the dispersion relation; whereas in specific examples one can derive the form of the power 
spectrum and precisely, which in the case of small speed of sound models turns out to be 
a factor of from saturating this bound. In any case, we recover the condition |^ < 10^ for 
a broader class of models, as anticipated. 

5 Prospects for observations 

We now discuss the observational constraints on the models discussed in previous sections. 
5.1 Non-adiabatic evolution 

The signatures of resonance models have been discussed in the literature, so we will simply 
summarize the relevant results. 

The strong-time dependence that allows us to violate the consistency condition also tends 
to produce large violations of scale invariance in the power spectrum. This is not surprising, 
as scale invariant correlation functions arise when each mode follows the same history. The 
time dependence of the couplings means that modes that cross the horizon at different times 
see different values of the couplings. It is this violation of scale invariance in the power 
spectrum that is the dominant signature of these models [32, 45] rather than the violation of 
the consistency condition for some intermediate regime of scales. This is not special to the 
specific realization of the model but is a generic feature of resonant models [45] . The reason is 
fairly intuitive. Once the approximate shift symmetry is broken, generically nothing forbids a 
contribution to the power spectrum at the same order as in the higher correlation functions. 
Then, since the spectrum is approximately Gaussian, the power spectrum will dominate the 
signal-to-noise ratio. 

One can try to avoid the above conclusion by introducing additional symmetries that 
forbid the large contributions to the power spectrum. Models of this type in which the 
bispectrum did contain most of the signal-to- noise were constructed in [46]. However, the 
models found there did not lead to a non-trival squeezed limit. The reasons are similar 
to what one finds for excited states with small c^, namely that the large interaction terms 
that increase the relative significance of the bispectrum are accompanied by derivatives that 
suppress the squeezed limit, even inside the horizon. 

For the case of dissipation, as we mentioned earlier, the scaling of the three-point function 
is not modified in the squeezed limit [13], however, one finds a peak at the folded configura- 
tion [12]. This is connected with the fact that particles are being produced during inflation 
[33-35]. This enhancement for the folded shape of the bispectrum occurs for any excited 
state, motivating the exploration of this (and other) shapes in the data. 

In summary, we are not aware of any model that produces a non-trivial squeezed limit 
and at the same time does not lead to a large signal in the power spectrum. 
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5.2 Bogoliubov initial states 

The excited initial states we have discussed were related to the Bunch-Davies state by a 
Bogoliubov transformation. As we mentioned earlier, the regime ^ 1 is ruled out by 
observations of the power spectrum. To see this recall that the power spectrum takes the 
form 

k'Uj:,m' = ^ (l + m\' + ^^/W + \^'^os{ej:)) . (5.1) 
As we explained, negligible backreaction on the background evolution enforces 

ks/kL < A-'/' ~ 10-2, 

and we can only significantly excite modes over two orders of magnitude. For k > k^, the 
Bogoliubov coefficients must decrease at least like /3g oc fc"" with n > 4. For generic values 
of 9^, this implies that the power spectrum scales like oc fc"^-" for k > k^, leading to a 
step in the power spectrum of primordial scalar perturbations. Observations of the cosmic 
microwave background alone by WMAP [1], SPT [2] and ACT [3], have reliably measured 
the angular power spectrum of temperature anisotropies over a range that is large enough to 
rule this out. We thus conclude that ^ 1 is incompatible with existing data. 

There appears to be a loophole in the above analysis for the case of Bogoliubov states 
when cos(6'g) = — 1 and /3g ^ 1. The power spectrum becomes independent of (3^ up to 
/3^2 corrections. However, to make use of this to avoid our conclusions, one would still need 
/3g 3> 1 for all observable modes. For k > k^ we have oc k~^ with n > 4. Therefore, 
to accomodate one more decade of highly excited modes to be compatible with observations, 
one needs /3g ^ 10^. However, this large occupation number is only compatible with the 
bound on the energy density of the excited modes, which takes the form ^ 0^ (see 

(3.21)), if k^^/ki ^ 10, in contradiction with the observation of the power spectrum. One may 
wonder whether other types of cancellations could occur for more complicated excited states. 
Again, this cancellation is most likely to occur for large occupation numbers, in which case 
our previous reasoning applies in the same manner. (See [22] for more details.) 

For the case of /3g "C 1, the constraints are weaker for two reasons. First, the contri- 
bution to the power spectrum is sub-leading relative to the scale invariant contribution, as it 
is nearly in the vacuum state. Second, the lower occupation numbers weaken the constraint 
on ks, allowing for the possibility of all the observed modes in an excited state with low 
occupation number ~ ^1. Nevertheless, the question remains if states with <^ 1 
could produce a measurable effect in the squeezed hmit. We now discuss two (in principle) 
possible routes for observation. 

5.2.1 Observing the squeezed limit in the CMB 

For the CMB, we can gain fairly reliable intuition for the signal-to-noise ration S/N by 
working with the correlation functions directly. In appendix B we derive the signal-to-noise 
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ratio for a direct measurement of the bispectrum, 



y ~2e/3,.^^xA,V:A^, (5.2) 

where A/'pix is the number of modes measured. A direct measurement of the bispectrum 
of primordial curvature perturbations is of course unreahstic, but incorporating effects of 
transfer functions, sky cuts, etc. necessarily decreases the signal-to-noise ratio so that this 
surely provides an upper bound. On the other hand, from the bound in equation (3.24) we 
notice that, for small occupation number, at best ks/k^ ~ 10^/3T^^^ so that we do not gain 
any signal by decreasing the occupation numbers to increase ks/k^. Therefore, the largest 
signal is expected^ for /33 < 10^^ and ki/ks ~ 10"^. Similarly, slow-roll inflation requires 
e < 10~^. Hence, using the number of observable modes in the CMB, A/'pix ~ ^max ~ (2000)^ 
and Aj^ = 2.2 x 10~^, one finds 

(I) r 

So even with the most optimistic choices of parameters, we see that no signal is expected with 
Planck [21]. 

We might hope to increase the signal at small by introducing interactions with large 
coefficients. We study this possibility in Appendix B for small Cg models. In the absence 
of an excited state, these models have an enhanced j^*?^^'''^*'^''^^ and can therefore lead 

to a detectable signature in equilateral configurations {ki ~ ^2 ~ k^). However, even in the 
presence of an excited state, there is no significant improvement for the signal-to-noise in the 
squeezed limit. Therefore, even allowing for stronger cubic interactions, these models cannot 
produce measurable deviations from the consistency condition in the CMB. 

5.2.2 Observing the squeezed limit in LSS 

Large scale structure surveys [4] have emerged as a powerful probe of squeezed limits of non- 
Gaussian correlation functions [49]. However, unlike the CMB, it is less clear what range 
of momenta contribute to the signal-to-noise for a given experiment and shape. For this 
reason, we will explore a forecast for a Euclid-like LSS survey [50], rather than using the 
simplified analysis we followed for the CMB. A detailed forcast is certainly beyond the scope 
of this paper and not even currently possible because the specifications of the mission are still 
evolving. However, using the specifications from the Red Book [51] as a guide for a simple 
forcast should allow us to draw robust conclusions especially because incorporating details 
would only decrease our idealistic signal-to-noise estimates. 

We expect the largest signal-to-noise ratio for a detection of the Bogoliubov initial states 
with low occupation numbers to arise from the the galaxy bispectrum. However, since we do 
not expect significant changes to our conclusions, we will focus instead on the measurement 

^We take this value as an upper limit consistent with the approximate scale invariance of the power 
spectrum. 
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of the 3D power spectrum 

Vg{k) = blPUk) . (5.4) 

Here Vm{k) is the underlying matter power spectrum and bg is the galaxy bias. 
The power spectrum is usually measured from the spectroscopic part of the survey. Euclid is 
expected to use a slitless spectrometer with a wavelength range 1100 nm < A < 2000 nm. The 
redshift measurement will rely on emission lines of galaxies, predominantly Ha, so that this 
wavelength range translates to a redshift range for the observed galaxies of 0.7 < z < 2.1. 
With a limiting flux and success rate of the spectrometer as given in [51], one expects a sample 
of 52.5 X 10^ galaxies over the survey area of the "required" 15, 000 square degrees. We will 
be optimistic and use the "goal" of 20, 000 square degrees rather than the "required" 15, 000 
square degrees for our forecast. 

Primordial non-Gaussianity has several effects on the galaxy power spectrum. The dom- 
inant effect is to introduce a scale-dependence into the bias 

bg{k, Jnl) = bgifNL = 0) + Sbgif^L) + Abg{k, f^L) ■ (5.5) 
According to the halo model, the galaxy bias is related to the underlying halo-bias by 

oo oo 

bg{kJ^L)= j dMbh{k,hL,M f^^^^'^^^ With ng= j dMn{M)Ng{M), (5.6) 

-^'^m in -^'-^m i n 

where Ng{M) is the number of galaxies expected in a halo with mass M, and n{M) is the 
number density of halos with mass between M and M + dM. The main lesson we take 
from these equations is that the scale-dependent contribution of the galaxy bias inherits its 
properties from the scale-dependent halo bias, which can be calculated as in [52-56] 



Mh{kjNL,M) 



(&i - 1)4 dHk.M) 
MM{k) [ D{z) ' dlogaM 



(5.7) 



where bi is the effective (Eulerian) Gaussian halo bias, 6c is the collapse threshold, and aM is 

1 /3 

the matter power spectrum smoothed on the scale R = {2MG /VLmH^) . Furthermore, D{z) 
is the linear growth function normalized to unity at redshift zero, M.M{k) = k'^T{k)W{kR) 
where T{k) is the transfer function in the conventions used by GAME, and W{x) = 3(sinx — 
X cos x) I x^ is the Fourier transform of a top-hat function. Finally, 

nk.M) = ^^2^^^ I (0jA^m(p)A^a/(I^ + p1)(CpC.<-,^_,^)'. (5.8) 

The resulting halo bias generated by the excited states discussed in Section 3.2 was computed 
in [41, 42]. We show it in figure 2 for a halo mass of 10^'^h~^Mso\, where the green and red 
lines represent the halo bias with and without our bound (3.24) enforced, respectively. 
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Figure 2: The green and red curves show the scale dependent halo bias as a function 
of k with and without the bound (3.24) enforced, respectively. We have chosen 
M = 10^2/^-1 Msoi, = 0.3, h = l Mpc-^, and slow-roll parameters e = r] = IQ-^. 



Let us now estimate what we can observe from a measurement of the power spectrum 
in our Euclid-like survey. We will assume a Gaussian likelihood function in the various 
momentum and redshift bins and calculate 

where Vg denotes a fiducial galaxy power spectrum. For the last equal sign, we have traded 
the sum over momenta for an integral, and we have used that the cosmic variance for the 
power spectrum measurement, up to a factor of 2, is just given by the effective number of 
modes that enter into the measurement (see e.g. [57]) 

^'ik,z^,) ^ 2(2vr)3 
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Figure 3: This plot shows Ax^ as defined in (5.11) as a function of with the 
phase and fc^ chosen to maximize the signaL The red and green curves show the 
results without and with the bound (3.24), respectively. The shaded region indicates 
the region that is ruled out by the power spectrum when (3.24) is enforced. 



Here Kff(^; Zi) is the effective volume. Especially for a fiux limited survey like Euclid it may be 
significantly smaller than the geometric volume [58].^ We will nevertheless use the geometric 
volume for our estimates, which will surely lead to an optimistic signal-to-noise estimate. We 
fix all cosmological parameters characterizing the background to their best-fit values [1], and 
use a fiducial power spectrum derived assuming a Bunch-Davies state. Keeping only the effect 
of the excited state on the scale-dependent bias, equation (5.9) then simplifies to 



(2 



k'dk 



1 + 



Abg{k,Zi 

hg{Zi) 



(5.11) 



In a more conservative treatment, one should treat the Gaussian bias as nuisance parameter. 
Instead, we model it using a Gaussian halo bias appropriate for the Sheth-Tormen halo mass 



^As shown in [58], the effective volume for different semi- analytic models of galaxy formation for Ha 
selection may be as small as a few per cent of the geometric volume for realistic Euclid flux limits and success 
rates. 
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function [59-61] and compute the effective galaxy bias from (5.6). We use tlie Sheth-Tormen 
lialo mass function for n{M), and use a simple fitting form for the expected number of galaxies 
in a halo with mass M [62] 



with iVgBiue = 0.7, Mq = 4 X lO^^Mgoi/i^^ and as = or as = 0.8 for masses below or 
above Mbiue, and we set M^^^ = lO^^Msoih'^ in equation (5.6). This roughly reproduces the 
results found for the effective Gaussian bias of an Ha sample typical for Euclid found for 
semi-analytic galaxy formation models in [58]. 

We split the survey into 12 redshift bins and set the limits of integration in (5.11) to 
^min.i = '^T^/V{z,i)^^'^ as wcU as fcmax.i = O.l/i Mpc~^ . We excite modes up to fc^ over a range as 
large as allowed by our bound (3.24) choosing as well as the phase to maximize (5.11), 
and set e = ?7 = 10"^. The result is shown in figure 3. 

We see that even though we have made a number of optimistic assumptions, the signal 
is undetectably small at least in the regime that is not already ruled out by power spectrum 
constraints. 

6 Discussion 

In this paper we have studied single-field inflation with excited states arising either from 
non-adiabatic evolution or as remnants from a pre-infiationary epoch. We have shown that 
all existing models in which the state is continuously excited during inflation can only vio- 
late the consistency condition for a finite range of momenta bounded by ks/ki < 10^. This 
implies that the violation of the consistency condition is unobservable in these scenarios. 
For resonant models and generalizations thereof the violation of the consistency condition 
comes with periodic features in the power spectrum, which provide the dominant signature in 
these models. Strongly dissipative systems give rise to a scale-invariant power spectrum, but 
the consistency condition applies roughly in the same fashion as in single-field models with 
Bunch-Davies states, which in practice means that the squeezed limit is not useful to distin- 
guish between the two and other corners of the bispectrum need to be explored, in particular 
folded configurations [12, 13]. For Bogoliubov initial states with significant occupation num- 
bers, but an energy density that is small enough so that the slow-roll conditions are satisfied 
(in particular S = <^ 1), we have also shown that the consistency condition can only be 
violated for a finite range of momenta bounded hj ks/ki ^ 10^. Thus in this case the violation 
of the consistency condition comes with a step in the power spectrum in contradiction with 
observations. For Bogoliubov initial states with small occupation numbers, the constraint is 
weaker. However, the three-point function is generated by gravitational interactions for these 
states. As a consequence, as we show in the form of signal-to-noise estimates, the violation 




(5.12) 
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of the consistency condition is too small to be observed in Planck or Euclid. 

What prevented us from extending our conclusions to models in which the energy density 
in the excited state, p^r, is larger than the kinetic energy in the inflaton (so that the condition 
5 <^ 1 is violated) is the potential breakdown of a controlled theoretical framework. In this 
case the bound ks/k^ < 10^ generically no longer applies, and a weaker bound on ks/ki might 
lead to observable consequences in the bispectrum without jeopardizing the scale invariance of 
the power spectrum over the observed range of scales. This regime therefore deserves further 
scrutiny, and one has to carefully study all sources of density perturbations and relevant scales 
in the n-point correlation functions, while ensuring that perturbation theory remains under 
control. 

Even though we may consider models when p^^ exceeds the kinetic energy in the infla- 
ton, we should always impose the so-called hack-reaction constraint, p,r < M^H^, since its 
violation would spoil inflation at least for some period of time. (In fact, using the null energy 
condition we can show more generally: jp^ < M^\H\7) What is no longer straightforward in 
the regime p-j^ ^ 0^ is how to relate this bound to the power spectrum. However, let us make 
two important observations about this regime. Firstly, if all the observed modes exit after p^^ 
has redshifted away (dropped below 0^), our conclusions remain unaltered. Secondly, since 
during the non slow-roll phase e ~ a"'*, we have ki/ko ~ (eo/^i)^''^ for modes with co-moving 
momenta A;o{i) = clo{i)H. This means that, for each decade of modes that exit during this 
phase the value of e in the subsequent slow-roll epoch decreases by a factor of 10^. Recall that 
at most two decades of excited modes with large occupation number can exit during the slow 
roll phase. Therefore, in order to accommodate the observed 10^ number of modes we would 
need e,^ < 10"'^ (taking eo ^ 1). This weakens the size of the (gravitational) interactions that 
produce the bispectrum for Bogoliubov states. Conversely, large field models with > 1 
cannot accommodate all three decades of (scale invariant) observed modes. We will discuss 
this and other aspects in detail in a forthcoming publication [22]. 

NOTE: While this work was being completed we were informed by A. Aravind, D. Lorsh- 
bough, and S. Paban of their related work [66]. 
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A A useful inequality 



In the main text, we were able to bound the ratio of scales ks/ki in terms of M^i\H\. In 
order to achieve an upper bound we must then find an upper bound on M^^\H\ in terms of 
H and fixed numerical coefficients. We will show this is possible in this appendix. 

Let us make the simplifying assumption that the kinetic term for vr is of the form A'^n'^. 
One can show that our result generalizes to arbitrary kinetic terms, but this will be unnec- 
essary for our purposes. In the EFT of inflation this kinetic term arises from two separate 
operators in the Lagrangian 

/bf4 

Ainetic = MlHd^it + 7r)9'^(t + vr) + -^[d^it + Ti)d^{t + vr) + 1]^ . (A.l) 

Only the first term in (A.l) gives rise to the gradient piece a~^(9i7r)^. As a result, one 
requires that M| > [18], otherwise we could send super-luminal signals using tt.^ Using this 
constraint we find 

A^ = 2M^ + Ml^\H\ > Ml^\H\ . (A.2) 

In order to translate (A.2) into a useful upper bound on Mp[|iJ| we must re-write A"^ in 
terms of known quantities. Because of the form of the kinetic terms in the action we know tt 

3 _ 1 

scales as A~^/c2a; 2. We also know the power spectrum is given by 

^ f^,c-/{oo ~ if) ~ 10-9 , (A.3) 

where we used ( = —Hn and that freeze-out occurs at u ^ H. Here we have assumed there 
are no large dimensionless parameters fioating around. This is clearly not true for the case 
of excited states with large occupation number, i.e. f3j: 1, which we cover below. We can 
now show the following inequality applies: 

< lA-%^u ~ H) . (A.4) 

Now let us return to the case of excited states considered in section 4.2 by modifying 
(A.3) to include non-trivial occupation numbers, i.e. 

.2 ^^.-3/ 



- ^4^p (^ -H)[l + 2|/3gJ^ + 2^|/3,^J2 + |/3,^J4cos(%)J , (A.5) 
where we assumed does not vary significantly with k. Therefore, the bound becomes 

^Pil^l < lAfc-'iu r^H)\l + 2\M' + 2,^iJT^cos(%)l . (A.6) 



^4 ^ 4^C P 



^Using causality and dispersion relation arguments, one can show M| > in flat space [63, 64]. 
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B Estimating signal-to-noise for the CMB 



In this appendix, we will estimate the signal-to-noise for CMB measurement that are expected 
from an excited state (in slow-roll). We will make this estimate by assuming we measure the 
curvature perturbations directly, rather than temperature fluctuations. This will serve as an 
upper bound on the signal-to-noise in the CMB. 

B.l Signal-to noise estimate in slow-roll models 

If we were to measure the bispectrum of ( directly, the signal-to-noise of this measurement 
would be given by 

In section 5, we showed that excited states with /3g > 1 are ruled out by observations of the 
power spectrum. Hence, we will work with Bogoliubov states in the limit of small and real 
such that [39] 

]_ _ ^i(ki+k2-k-i,) / ki, _ g«(fci-fc2+A;3)/fc* _ ^i{kz+k2-ki)/ki, 

^^^^ ki + k2-h ^ ki-k2 + h k3 + k2-ki ^ ' ^'^'^^ 

(Here we have dropped the 0{(3'i ) contribution because it is unobservable.) 

The dominant contribution to the signal-to-noise comes from the flattened limit ^2 = 
^3 ~ ki (and equivalent configurations). We may therefore simplify the bispectrum in this 
folded limit to get 

1 f76 1 1 _ pi(k3+k2-ki)/ki. 

Evaluating the signal-to-noise analytically is possible, but the full expression is not very 
insightful. We are most interested in its parametric dependence, which we can understand 
with the aid of some approximations. First of all, we notice that the exact folded limit is cut 
off by the oscillatory term when X3 ~ 1 — Xi -|- x^, where x^, = k^/k^,. Using a hard cutoff we 
can then evaluate the integral over X3, while the X2 integral receives its largest contribution 
form X2 ~ 1. Therefore we have, 

f^i f d% 3! 1 

dx2j- — (B.5) 



NJ 4 MpV^Asy (27r)3 (47r2) 7 ^ {1 - X2)^x 
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where A/J^ix is the number of independent modes we can measure (i.e the volume of the ki 
integral), and we have used ~ a^max ~ k^/ks to measure the maximum ratio of momenta 
where we have the largest contribution to the squeezed limit. 

B.2 Signal-to-noise estimate in small-c^ models 

In small Cg models there is potential increase in the bispectrum due to larger interactions. 
Specifically we will consider 

i/,«='^/<iW*^, (B.7) 

cj J 

in the limit Cg ^ 1, that follows from the EFT of inflation [18]. Although the strength of 
the interaction is enhanced by a factor of 1/c^, we also expect these will be suppressed in the 
squeezed limit due to type of derivatively coupled interactions. 

We consider again Bogoliubov states with <^ 1 due to the power spectrum constraint. 
Following the same steps as for the slow-roll case, the dominant contribution to the signal-to- 
noise comes from the folded limit k2 = k^ — ki. We may, as before, simplify the bispectrum 
in this limit and get 

^ A^m - kl - kl) + permutations) ^^^ 1 - e'{k,+k,-k,)/k 
B{k,, k„ ks) - ^3^P3 Re[ ] , (B.8) 

where we used that, <^ 1, the mode functions for C, take the form 

Cs = —■ 7^7 (1 - iCskT)e''''^^ . (B.9) 

Again we wish to calculate the signal-to-noise: 
SV f^l .2 f d% 3! 



and we find 



Jl/2 Jl-X2+Xi, ^^^^ 



nJ 16cf ^ J (27r)3 (47r2 



1/2 



Notice we do see an enhancing factor of Xj,"^, however it is easy to see that this integral 
does not get contributions from X2 ~ 1 (or equivalently X3 — 0). As a result, there is no 
improvement for the signal-to-noise in the squeezed limit for small speed of sound models. 
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